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1 .  INTRODUCTION 

There  are  three  versions  of  the  finite  element  method:  the  In¬ 
version,  the  p-version  and  the  h-p  version.  The  h-version  is  the  standard 
one,  where  the  degree  p  of  the  elements  is  fixed,  usually  on  low  level, 
typically  p  =  1,2,3  and  the  accuracy  is  achieved  by  properly  refining  the 
mesh.  The  p-version,  in  contrast,  fixes  the  mesh  and  achieves  the  accuracy 
by  increasing  the  degree  p  of  the  elements  uniformly  or  selectively.  The 
h-p  version  is  the  combination  of  both. 

The  standard  h-version  has  been  thoroughly  investigated  theoretically 
(see  eg.  [1],  [9],  [19]  and  others)  and  many  program  codes  are  available, 
both  commercial  and  research.  The  p-version  and  the  h-p  version  are  new 
developments.  There  is  only  one  commercial  code,  the  system  PROBE  (Noetic 
Technologies,  St.  Louis).1  The  theoretical  aspects  have  been  studied  only 
recently.  The  first  theoretical  paper  appeared  in  1981  (see  [6]).  See 
also  [2],  [5],  [7],  [10],  [11],  [14]  for  most  recent  results.  For  the 
numerical,  computational,  impleraentational  and  engineering  aspects  of  the 
h-p  version  we  refer  to  [3],  [20],  [21],  [22]. 

The  classical  form  of  the  error  estimate  for  the  h-version  with 
quasiuniform  mesh  is 

(1.1a)  |u  -u  |  S  C(p)hn'1|u| 

U  *“  H  ( 0)  U  fTU) 

where 

(1.1b)  n  =  min(k,p+1 ) 


1  In  addition  there  is  code  FIESTA  for  solving  3  dimensional  elasticity 
problems  having  p-version  features  but  using  only  IS  p  S  4. 


V-V- 


«•_  V. 


*  M  • ' 


'  y  j  vlv  1  vlvlvlfoy]  .-rl  vlv? 


•jVi 


and  the  constant  C(p)  depends  on  p  in  an  unspecified  way.  See  eg. 
[1],  [9],  [19]  and  others. 


The  main  purpose  of  this  paper  is  to  analyze  the  h-p  version  with  a 
quasiuniform  mesh  and  uniform  p  and  get  an  error  estimate  which  is 
simultaneously  optimal  in  both  p  and  h.  We  show  that  the  estimate  (1.1) 
can  be  written  in  the  form 


(1.2) 


lu0_UFE^  1 

U  H  ( fl) 


hn-i 

p  H  (Q) 


with 


n  =  min(k,p+1 ) 


and  C  independent  of  h,p  and  u.  We  will  also  prove  estimates  for  the 
h-p  version  when  the  solution  has  singularities  in  the  corners  of  the 
domain  and  in  the  case  when  essential  (Dirichlet)  conditions  are  prescribed 
but  are  not  in  the  subspace  of  finite  elements.  Finally,  we  will  present  a 
numerical  example  illustrating  the  applicability  of  the  developed  (asymp¬ 
totic)  theory  in  a  range  of  h  and  p  used  in  practice. 


2. 


THE  NOTATION 


For  8  c  r2  a  polygonal  domain,  x  =  (x1,x2)  €  R2,  we  let  L2(8)  = 
H°(8),  Hk(8),  Hq(8),  k  >  0  integer,  denote  the  usual  Sobolev 
spaces.  For  u  €  Hk(8)  we  denote  by  |u|k  ^  and  the  usual  norm 

and  seminorm,  respectively.  For  k  £  0  nonintegral,  we  define  Hk(8) 

and  |'|  by  the  K-method  of  the  interpolation  theory  [8],  If  I  is 

K  )  U 

an  interval  or  a  segment,  then  we  define  Hk(I),  |*|,  _,  k  £  0 

analogously. 

Given  p  >  0,  let  R(p)  =  {(x1tx2)  |  Ix^  <  p,  |x2|  <  p}.  For  any 
8  c  r2  we  win  denote  pn  =  sup{diam(B)  |  B  a  ball  in  8}. 

The  set  of  all  algebraic  polynomials  of  degree  (total)  less  than  or 
equal  to  p  on  8  will  be  denoted  by  Pp(8).  By  P2(8)  we  will  denote 
the  set  of  all  polynomials  of  degree  less  than  or  equal  to  p  in  each 
variable  on  8.  For  r  c  R2  a  straight  segment,  we  define  Pp(D  as  the 
set  of  polynomials  on  r  of  degree  less  than  or  equal  to  p  in  3  (s 
being  the  length  parameter  of  D. 

Let  <  >  0.  Then  by  HpER(R(te))  c  Hk(R(<))  we  denote  the  space  of 

1  p 

all  periodic  functions  with  period  2k.  By  Tp(R(<))  (Tp(R(ic))  we  denote 
the  3pace  of  all  trigonometric  polynomials  of  (total)  degree  (degree  in 
every  variable)  less  than  or  equal  to  p. 


'■'uCe.  io  For 

~  — t 

^TiS  CRA&I 
or ,C  TAB 

u  "’r.rc  :•  cf-d 

— 1 

_ 

j 

By 

Di-t  ib;.tio:  /  j 

Awdiiaiv-ry  c 

'•c.'es  i 

The  internal  angle  at  is  denoted  by  We  allow  the  possibility 

that  -  ir  or  2ir.  The  case  ui^  =  2ir  describes  the  slit  (cracked)  domain 
while  the  case  =  it  is  introduced  to  deal  with  the  abrupt  change  of  the 
type  of  the  boundary  condition  or  with  nonsraoothness  of  g  or  b  at  the 
corresponding  vertex.  When  fl  is  stated  to  be  a  Lipschitz  polygonal 
domain,  then  it  will  be  assumed  that  <  2ir,  i  =  1,2,...,n. 

Let  Hg(SJ)  =  {v  €  v  =  0  on  r1  } .  For  u,v€H1(Q)  we  let 

(u,v)q  q  =  j  uvdx,  (u,v)1  n  *  /  (Vu  •  Vv  +  uv)dx.  We  interpret  now 
(3.1)  and  (3.2)  in  the  standard  variational  sense  namely  we  seek  u  € 

H1 ( Q)  so  that 

(3.3a)  u  =  g  on  r1 


and 


(3.3b) 


(u’V)i,a  =  (f*v)o,p  +  1 


bvds 


holds  for  all  v  6  Hq(£2). 

We  will  assume  that  the  solution  u  of  (3.1)  and  (3.2)  is 


(3.4) 

where 

(3.4a) 


(3.4b) 


(3.4c) 


U  =  U^  +  U£  +  Ug 


u,  e  h  ](n)  n  h1(q). 


Up  €  H  *(Q), 


r.1 


'3  ’  J,  aiu3,l  f 


(n>. 


u3,  i  =  ri 1  I  log  r,- 


k1  >  1 


kp  >  3/2 


(3-4d) 


Y. 

1 


(8. )x, (r.  ) 


where  r^,  9^  are  polar  coordinates  with  respect  to  the  origin  located  at 

the  vertex  A^,  >  0,  >  0  integer,  «P i ( 0 i )  is  an  analytic  function 

in  and  x ^ ^ r ^ )  is  the  C*  cut-off  function  so  that  u^  ^  =  0  for 

ri  >  Pj  >  0,  pi  sufficiently  small. 

The  form  (3.4)  is  the  typical  form  of  the  solution  of  (3.1)  (3.2)  (and 
of  a  system  of  second  order)  (see  eg.  [4],  [12],  [16]).  The  assumption 
that  k  >  3/2  is  usually  satisfied  in  practice  and  hence  is  not  a  severe 
restriction. 


3.2.  The  finite  element  method 

Let  M  =  {Th ),  h  >  0  be  a  family  of  meshes  T 11  =  {S^}  where  c 

Q  is  an  open  triangle  or  parallelogram.  Let  h  .  =  diam(Sh)  and  p  h  be 

Sn  S 

as  defined  in  Section  2.  We  shall  assume  that  the  family  {Th}  is  regular 
in  the  sense  that  there  exist  positive  constants  o,  t  independent  of  h 
such  that  for  all  sj  6  Th,  Th  €  M 


(3.5a) 


max 


h 


(3.5b) 


(3.5c) 


T 


0. 


(Condition  (3.5b)  is  obviously  the  condition  of  quasiuniformity  of  the 

mesh).  Further  we  assume  that  with  Th  =  (S^},i  =  1,2,...^^, 

m, 

—  ^  — h  —  — w 

3  =  J  S  and  that  each  pair  S” ,  S'),  i  /  j  has  either  an  entire  side 

i  =  i  1  J 


or  a  vertex  in  common,  or  has  empty  intersection. 
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Let  F*?  be  an  affine  mapping  with  Jacobian  having  positive 

<J 

determinant  which  maps  Sj  onto  the  standard  square  Q  =  (-1,1)  *  (-1,1) 

when  Sj  is  a  parallelogram  and  onto  the  standard  triangle  T  =  {(x.,^)  | 

-1  <  x1  <  1 ,  1  <  *2  <  xi  ^  when  Sj  is  a  triangle.  Let  now  l/p(D)  c 

hVj))  denote  the  set  of  functions  u  such  that  if  u  ,  denotes  the 

on 

5i 

restriction  of  u  to  S?  €  then  u  ^  o  (Fb)_1  €  Pp(Q)  if  is  a 

Si 

parallelogram  and  u  h  o(Fb)-1  €  Pp(Q)  if  is  a  triangle.  We  will 

Si 

then  write  u  y,  €  Pn(S^)  and  u  €  l/b(Q).  Furthermore,  we  let  l/b(ft)  = 

5V  P  1  P  P 

l/p (52)  n  hJj(Q). 

The  mesh  Tn  on  52  induces  a  partition  =  {^i  j),  j  = 

1,2,...m(i)  of  rt,  i  =  1,...,n  .  Denote  by  Nb  , ,  j  =  0, 1  , . . .  ,m(  i )  the 

nodal  points  of  L^1  (i.e.  the  end  points  of  Y1?  ^).  We  let  l/!?(r,)  c 

hVf^)  be  the  set  of  functions  u  such  that  the  restriction  u  ^  of 

Yi.j 

Vm  o  I-,  V* 

u  on  Y^j  is  a  polynomial  of  degree  £  p.  Moreover,  (  ri  ^  c  ^  ri  ^ 
will  denote  those  polynomials  that  vanish  on  N1?  . ,  j  =  0, 1 , . . .  ,m(  i ) . 

1  fJ 

Let  gb  €  U  ^ d ^ f i ^  be  tbe  aPProximati°n  of  g  (see  (3.2))  described 
p  r.cr  p 
below.  1 

The  h-p  version  of  the  finite  element  method  consists  now  (for  given 
p  and  h)  of  finding  uj  €  Vp(fl)  SUCh  that 


(3.6a) 


(3.6b) 


=  g*  on  r1 


(u,v)1fJJ  =  (f,v)0ffl  +  j  bvds 


holds  for  all  v  €  l7p(ft) 


To  define  g£  we  denote  by  gr  the  restriction  of  g 

f,  c  r1  and  assume  that  gr  €  Hr(r,.  )  with  r  >  1. 
l  li  1 


'.'.**** 


/  aV-  V  A A  .  -  „  •  A  A  A  a  A  -  A  /vV  /V*  - .  w  v  v  *.  -  -  *.« 

*.  ’  •  .  i  >  _  ■  'ji  *>  V  .'j*  'j  *  »  *  *  •  ‘  .*  >  A  V  *  /  '  * 


(3.7c)  |  (g„  r  )’w'ds  =  f  g'w'ds,  r,  c  r1 

r .  “  * 1  i  r  1 

i  x  x 

holds  for  all  w$  ^ p ( T ^ ) - 

Remark.  If  we  restrict  (3*7b)  to  j  =  0,m(i)  only  ( N*-1  n  =  A., 

Nj  m(i)  =  A^+^),  then  (3.7b)  is  satisfied  as  a  consequence  of  (3.7c). 


4.  THE  CONVERGENCE  OF  THE  h-p  VERSION:  THE  CASE  OF  THE  SOLUTION  u€HK(Q) 
In  this  section  we  will  analyze  the  rate  of  convergence  of  the  h-p 
version  when  the  solution  of  (3.1),  (3.2)  has  the  form  (3.4)  with  =  0. 

4.1.  Basic  approximation  results 

We  present  here  some  approximation  results  which  will  play  an 
essential  role  later. 

Lemma  4.1.  Let  S  =  Q  or  S  =  T  be  the  standard  square  or 

triangle.  Then  there  exists  a  family  of  operators  { Trp ) ,  p  =  1,2,3,..., 

tt  :  Hk(S)  ■+  P_(S)  such  that  for  any  0  <  q  <,  k,  u  €  Hk(S) 

P  P 

( 4 . 1  a )  |u-ir  u|  s  <  Cp"0<-q)|u|k  s  ,  k  >  0 

(4.1b)  |(u-y)(x)|  *  Cp*(k_1)  lujk  s  ,  k  >  3/2,  x  6  S 

where  we  denote  Pn(S)  =  P^(S)  for  S  =  Q  and  P  (S)  =  P„(S)  for  3  = 

T.  The  constant  C  in  (4.1a)  (4.1b)  is  independent  of  u  and  p  but 
depends  on  k. 

Moreover,  if  u  (  Pp(S),  then  Trp(u)  =  u. 

Proof.  The  proof  of  this  lemma  is  an  adaptation  of  the  proof  given  in 
[5].  Hence  we  will  only  outline  the  proof. 

Let  rg  >  1  so  that  S  c  R(rQ) .  Since  S  is  a  Lipschitz  domain, 
there  exists  an  extension  operator  T  mapping  H^tSI  into  Hk(R(2rg)) 
such  that 

Tu  =  0  on  R(2r0)  -  R(J  rn) 


(4.2a) 


(4.2b) 


k,R(2r0) 


1  2 


9  Tu  I 


*  CK,S 


where  C  is  independent  of  u.  For  a  concrete  construction  of  T  we 
refer,  for  example,  to  [4],  [18]. 

Let  $  be  the  one-to-one  mapping  of  R(j)  onto  R(2rg): 


(4.3) 


R(2rQ)  5  x  =  (x1  ,x2)  =  <S>(  O 


=  (2rQ  sin  ^  ,  2rQ  sin 


with  (5, ,§2>  -  5  €  R(j). 
Further,  we  let 


(4.4) 


R  =  [R(|  rQ)]  cR(|) 


where  $  1  denotes  the  inverse  mapping  of  <J>.  Let  v  =  Tu  and 


(4.5) 


VU)  =  v(*(C)). 


Because  of  (4.2a)  we  easily  see  that 
(4.6)  Supp  V(£)  c  R. 

In  addition  it  can  be  readily  seen  that 


(4.7a) 

V  6  h£er(R(it)), 

(4.7b) 

|v|k,BU)  <  cMk,s  . 

(4.7c) 

vu) 

is  a  symmetric  function  with  respect  to  the 

lines 

5i  "  ±  f.  1  =  ^ 

Let  us  expand  the  f unction  V  in  terms  of  its  Fourier  series 


S.'S"'  \\V 

-A.*.*...' 
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VU..S?)  =  I  l  a  e 

J  —  2=-»  J 


For  any  p  >  1  we  define 


i)  for  S  =  Q: 


(4.8a) 


I  I  aio  e 

!jUp  |h|^p  j 


i(J51+U2) 


ii)  for  S  = 


(4.8b) 


a  .  „  e 


iCj^+Ug) 


|J  1*1  1Up  Jl 


Then  quite  similarly  as  in  [5]  we  have  for  0  £  q  £  k 


(4.9a) 


(4.9b) 


lv“*pVlq,R(w)  i  cP'0<',)|“lk,s 

|(V-*  V(5))|  S  Cp"(k"1)|u|.  _  , 


k  >  0 


k  >  3/2. 


Because  (TTpV)  (4>-1  (x) )  €  Pp(S)  and  *  is  a  regular  mapping  of  R(rQ)  (rQ 
<  -|)  on  S,  (4.9)  yields  immediately  the  lemma. 

Let  us  quote  now  the  following  scaling  result. 

Lemma  4.2.  Let  ft  and  fth  be  two  open  subsets  of  R°  such  that 
there  exists  an  affine  mapping  F(x)  =  B(x)  +  b  of  ft*1  onto  ft  and  F(fth) 
=  ft.  Let  diam(ft)  =1,  Pq  =  K,  diam(fth)  =  h,  =  Kh.  If  the 

function  v  €  Hm(ft),  m  >  0  integer,  then  v  =  v  o  F  €  Hm(ft!l)  and 


(4.10a) 


5  Ch 


(4.10b) 


S  Ch 


I  m,  fl 


M  nh 

m ,  fl 


where  C  depends  on  K  and  K  but  not  on  fl,  h,  v. 

For  the  proof  see  [9],  Theorem  3.1.2 

The  estimate  of  the  error  of  the  approximation  of  g  by  is  given 

in 


Lemma  4.3.  Let  r>1,  0  £  t  £  1 ,  p  >  1 ,  then 


(4.11a) 


u-«Ji  „n  i  c  Vt  igi  h 

p  r'Ti,j 


(4.11b) 


v-t 


lg-gp®t,r.  -  C  r-t  ®glr,r 


where 

(4.11c) 


v  =  min(r  ,p+1  ) 


and  C  is  independent  of  g,  p  and  h. 


The  proof  is  given  in  [13].  The  main  idea  is  to  expand  g'  in 

Legendre  polynomials  on  every  .  of  the  partitioning  of  r,  induced  by 

1  t  J  1 

the  mesh  T  ,  prove  (4.11)  for  r  and  t  integral  and  by  the  interpola¬ 
tion  argument  obtain  (4.1 i)  in  full  generality. 

Let  us  prove  now 

Lemma  4.4.  Let  31"1  and  S  be  the  triangle  or  parallelogram 
satisfying  the  conditions  of  Lemma  4.2.  Then  for  any  u  €  H'<(fl) 
corresponding  to  the  function  u  €  HK(fl^),  k  >  0  we  have 
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(4.12) 


.inf  |a-p|kQ  <  Ch^|u| 

P€Pp(Q)  k’°  k'Qh 


where  u  =  min(p+1,k)  and  C  depends  on  K,  R,  k  but  is  independent  of 
p  and  u . 


Proof.  For  k  =  0  the  result  follows  immediately  from  Lemma  4.2 
taking  p  =  0.  Hence  let  k  >  0.  Assume  fir3t  that  k  is  an  integer. 


.  inf  |G-p|  £  *  inf  l|G-p|  $  *  I  I ^ I t  q  +  l  IpI t  ^ 

p epp(fl)  k,n  p€pp(fl)  i-u+1  i  =  u+1 


where  £  =0  for  k  <  u  +  1.  Using  Theorem  3.1-1  of  [9],  we  see  that 

i  =  u+1 


.  inf  |u-p |  <  C  l  |u| 

peP  (0)  *  i-U 

r 


<  C  1  h1  1  |u|  (by  4.10b) 

i-u  i.n 


<  Ch^Jul 

k,n 


and  (4.12)  is  proven  for  k  integer.  For  general  k  we  use  an 
interpolation  argument. 

Let  us  prove  now 


Lemma  4.5.  Let  Sn  be  a  triangle  or  parallelogram  with  vertices 
satisfying  conditions  (3.5).  Let  u  e  Hk(Sl1).  Then  there  exists  a 
constant  C  depending  on  k,  t,  o  but  independent  of  u,  p  and  h  and 
a  sequence  zp  (  Pp(Sh),  p  =  1,2,...  (see  def.  of  Pp(Sh)  in  Section  3-2) 
such  that  for  any  0  i  q  £  k 


war®. 


(4.13a) 


h  *  C  M  h 

P  q,Sn  p  q  k ,Sn 


k  >  0 


(4.13b) 


|(u-zjj)(x)|  1  C  |u| 


k  >  3/2,  x  €  Sr 


(4.13c) 


u  *  min(p+1  ,k) , 


If  k  >  3/2,  then  we  can  assume  that  Zp(A^)  =  u(A^). 


Proof.  Let  ir  be  the  operator  introduced  in  Lemma  4.1.  Define  now 

itj:  Hk(Sh)  *  Pp(Sh) 


so  that 


TTpU  =  (lp(u  o  F*1))  o  F 


where  F  is  the  linear  mapping  of  Sn  onto  T,  respectively  Q  (see 
Section  3.2).  Denoting  u  =  u  o  F-1  we  get  from  Lemma  4.1  and  4.4  for 
q  ^  k 


(4.14) 


lG"V*q.S  =  3(u-p)  -  yu-p)Iq>s 


<  Cp"(k~q)  A  inf  |u-pi 

P€Pp(S) 

-(k-q)  u-1  ,  , 

<  Cp  h  H  h. 

k ,  3 


Combining  (4.14)  with  Lemma  4.2  we  get  for  0  £  m  Z  q  i  k 


i  h  |  s  rh^-mD~  ^ i  i 

lu-^m.s*  *  Ch  P  M,  h 


and  hence 


(4.15) 


lu'Vl  „h  -  chp'qp“u‘q)|u|  . 


Now  analogously  for  k  >  3/2  and  x  e  S 


(4.16) 


(u-ir  u)(x)|  *  Cp*(k_1)  ^  inf  |u-p |  _ 

P  P€PP (S) 


i  Cp'(k*1)hll‘1  Ju| 

k,S 

and  (4.13)  is  proven. 

If  k  >  3/2,  then  we  can  modify  Zp  by  a  linear  function  if  S  is 
triangle  or  a  bilinear  function  if  S  is  a  parallelogram  30  that  z£(A^)  = 
u(Ai).  Using  (4.13b)  it  can  be  readily  seen  that  (4.13a)  will  hold  once 
again  for  this  modified  function. 

The  proof  of  the  following  theorem  is  a  modified  version  of  Theorem 
4.1  in  [5]. 


Theorem  4.6.  Let  u  be  the  solution  of  (3. 1-3. 2),  u  €  Hk(Q), 
k  >  3/2  and  for  T1  c  T1  let  gj^  €  Hr(Ti),  r  >  k  -  V2,  where  gj  is  the 
restriction  of  g  to  T^.  Then  for  each  p  >  1  and  h  >  0,  there 
exists  <Pp  €  ^p(n)  such  that 


(4.17a) 


on  r1 


(4.17b) 


lu-tpp li ,q  ^  c  k-i  (Hu^k,n  +  £  ^g^r,r.) 


(4.17c) 


u  =  min(p+1 ,k) 


where  gq  is  defined  by  (3.7)  and  C  i3  independent  of  u,  p,  h,  and  Th. 


First  we  will  introduce 
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Lemma  4.7.  Let  S  =  Q  or  S  =  T  and  let  Y  =  A2  be  a  3ide  of 
S.  Let  ijj  €  Pp(Y)  such  that  (KA^)  =  0,  i  =  1,2.  Then  there  exists  an 
extension  v  €  pp(  S) ,  v  =  4)  on  Y,  v  =  0  on  3S-Y  and 


(4.18) 


I'll, 5  *  C|*|y 


where  the  constant  C  is  independent  of  p  and  \|i. 

The  lemma  is  a  special  version  of  Theorems  7.4  and  7.5  presented  in 
the  Appendix. 

Proof  of  Theorem  4.6.  Let  {s!?}  *  Th.  Then  by  Lemma  4.5  there 
exists  zpji  €  pp(S^)  such  that  zp  ^  =  u  at  all  vertices  of  S^.  Let 
now  Y11  =  Sj  fl  3^  and  let  N., ,  N2  be  the  end  points  of  Yh.  Then  zp  j  - 
zp,Z  =  wj,z  i3  a  polynomial  on  Yh  of  degree  at  most  p,  and  w*?  a( )  = 
0,  i  *  1,2.  We  now  map  S*?  (J  onto  U  Sj  by  a  continuous  linear 
mapping  F  where  Sj  and  are  congruent  images  of  Q  or  T,  suitably 

placed  as  shown  in  Fig.  4.1. 

Using  the  notation  used  in  the  proof  of  Lemma  4.5  we  get,  by  Lemma  4.5 


|u“zp  j  I,  s  ^  C  luL 

P'J  pK  1  k, Sj 

|u-z  oli  q  is  analogously  bounded.  Hence  on  Y  we  get  by  the  imbedding 
P  f  •  f 


theorem 


i-h  ■  .  « "  ''h  •  ~h  ■ 

|wjA.T  -  '“‘Vjll.Sj  *  iJ'2p,ii.,sf 


-  c  npr  (N,  j,  *  M,  h> 


v  V  'J  *  »  *>  mM  'V'J  •  h  *  *  O  n  *>  “ 


•  •  ►  •  »  •  a  <  «  / 

«•  ^ '  *»  •  »  "  »  ‘  w  a  *  4  «  *  •  • ■  *  «  “  »  •  *  a  w  la 


Fig.  4.1.  Scheme  for  the  map  of  two  neighboring  elements 


Applying  Lemma  4.7  there  exists  $£Pp(Sj)  so  that 


1*1,  s 

<  1 

Wj.i|l/2,Y 

A 

*  = 

Qi,i 

on  Y 

and 

A 

^  =  0 

on 

9Sj  -  Y. 

Hence  we  can 

modify  zi? 

P  ! 

J 

to  z5,j 

so 

that  zj}  ,  s 
P  t  J 

=  zp  ^  on  Y*1  and 

i  i 

Z?,J'  U  i  sh 
,Sj 

hu-l 
-  c  TT 

p 

(1 

^•sj 

J 

M 

k 

h>- 

,SH 

Repeating  this  process  we  construct  4  similarly  on  each  S*?. 

P  >  J  J 

Defining  so  that  its  restriction  on  Sj  is  2p(j  we  get  3p  £  l/p(£) 


i.* 

s* 

•  M 


|u- 


-h| 

<Pp  It 


hy_1 


fc-La- 


,  rm  . 


’-a  j*.  r.jp.s*,.  *«r.j 


jfcJ! 


*  o  “J* O’ * “ - • *  '  • *  .*«* 
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Finally  if  3S1?  fl  r1  =  Th  /  0,  we  have  to  modify  z j?  .  so  that 

J  P  t  J 

■zrl  ■  =  g1}  on  Yh.  Using  (4.11a)  and  realizing  that 
P  *  J  P 

1  Isl2  H  <  Ul2  r 

J  r.Y?  .  r,Ii 

i.J 

we  can  proceed  quite  analogously  as  before  and  complete  the  proof. 

Remark.  By  the  imbedding  theorem  we  have  |g  |  y  r  <  Iulk  Q  and 

*' 1  i  ' 

hence  the  second  term  in  (4.17b)  can  be  omitted. 


4.2.  The  approximation  results  for  1  <  k  <  3/2 

In  the  previous  section  we  analyzed  the  case  when  the  solution  u  of 
(3-1)  (3.2)  belongs  to  Hk(R),  k  >  3/2.  We  will  now  analyze  the  case 
when  u  €  Hk(a),  1  <  k  <  3/2  and  g  =  0.  In  addition,  we  will  assume 
that  fl  is  a  Lipschitz  domain. 

As  shown  in  [4],  given  any  t  >  0  and  k  >  1,  the  function  u  can  be 
decomposed  so  that 

(4.19)  u  =  vfc  + 

vfc  6  Hj(Q) 

?  Hk(n)  n  hJ(q) 

and  for  any  k  >  q  >  1 


(4.20b) 
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Let  2  >  k  >  3/2,  and  1  £  q  <  3/2.  Then  by  Theorem  4.6  there 
exists  <Pp  €  l/£(q)  such  that 

cpp  =  0  on  T1 


We  remark  that  the  assumption  that  Q  is  a  Lipschitz  domain  was  used  in  the 
proof  of  decomposition  (4.20).  (4.21)  shows  that  in  Theorem  4.6  we  can 

replace  the  restriction  k  >  3/2  by  k  >  1  provided  that  g  =  0.  In  fact, 
we  need  less  namely  that  g|p^  €  Hr(I\),  c  r1  ,  r  >  1. 
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4.3.  The  rate  of  convergence  of  the  h-p  version  of  the  finite  element 
method 


We  will  prove  now 


Theorem  4.8.  Let  u  €  Hk(ft),  k  >  1  be  the  solution  of  (3-1) 
(3.2).  Assume  further  that  g  is  such  that 


U  =  u-|  +  Up 


k1  —  i 

^  €  H  1  ( £2)  n  Hq( Q) 


up  €  H  ^(ft), 


k2  >  3/2 


and  that  ft  is  a  Lipschitz  domain  if  k1  £  3/2.  Let  u£  be  the  finite 

element  solution  of  (3.1 )— (3-2)  as  defined  in  Section  3-2,  then 


(4.22a) 


(4.22b) 


HS'i.a  *  c(k)  lulk, 


k  =  min(k1,k2) 


(4.22c) 


u  =  min(p+l  ,k) 


where  C  is  independent  of  u,  h,  p  but  depends  on  ft,  t,  a. 


Proof.  If  g  =  0  then  (4.22)  follows  immediately  from  Theorem  4.6 


and  (4.21 ) . 


If  g  4  0,  then  denote  by  (Jp  the  exact  solution  of  the  problem 


(3-1 )— (3-2)  when  replacing  g  by  gp*  Denoting 


U)  =  u  -  Up  we  see  that 


-  iiui  +  ii)  =  0 


'  *  •  * •  •  « « •» 


«  •  *  »  “  »  *  m  *  jt  *  i  ",  *  4*. 

Iv  .V.V 4*.  -V-V- 


„  *-  A  A  .  A.  -u  A  *, 


t  / 
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3u>  A  „2 

t—  =  0  on  r 

an 


u  =  g  -  gp  on  r1 


h  z 

By  Lemma  4.3  we  have  M  .  <  C  - 17  |u|  ,  where  v  =  min(r,p+1)  and  r 

l4.r  “  pr~/2  r ,  r1 

=  k  -  \  by  the  imbedding  theorem.  Because 

Mi, a  -  inrM,,a 


over  all  v  €  ^(8)  such  that  v  =  ui  on  F1  ,  we  have 


Mi, a  S  cl“l  j  i  c  —  K.tr 

4. "  P 


3y  Theorem  4.6  and  the  basic  properties  of  the  finite  element  method  we  get 
for  any  €  t/£(8). 


•  h  ,,h  1  .  „  1  h  tIh  I 

lu  -U  ,  „  <  C  m  -U  I,  , 

1  p  p 1 1  ,8  -  ,vp  p'l  ,1 


i  c(  M-^pli  ,£!  *  lu-upl,  ,Q> 


hu_1 

^  'TPTK.: 


and  Theorem  4.3  is  proven. 


4.4.  Optimality  of  the  asymptotic  rate  of  convergence 

In  this  section  we  will  prove  that  the  estimate  in  Theorem  4.3  is 
optimal.  To  do  so  we  will  use  the  concept  of  the  n-width.  For  details, 
see  eg.  [17].  Denote 


Dn(H]U),  Hk(Q)) 


i  m 


rk 


ini  y  Li  V  I 


Sch’cq)  u€Hk  v€Sn 

"  Ml  -"I 


1,0 


dim  S 


n 


'k,0 


the  n-width  in  the  sense  of  Kolmogrov.  Then  by  Theorem  2.5.1  and  2.5.2  of 
[1  ]  we  have 


:4.23) 


D-CH^n),  Hk(0))  >  Cn'1/2(k-1  ) . 


Let  us  now  compute  the  dimension  of  the  space  l/p(Q)  in  terms  of  p  and  h. 

The  number  of  elements  is  of  order  o(— j).  Over  each  element  we  have 

h  2 

0(p2)  polynomial  basis  functions.  Hence,  n  =  dim  l/p(0)  <,  C  Hence 

h 

for  p  +  1  >  k  we  have 


h  h  k_1  -(^> 

(U‘24)  lU~Up  ll,Q  <  C(k^o}  Mk.Q  ^  C(k)n  M 


■k,8 


k,fl* 


Comparing  (4.24)  with  (4.23)  we  see  that  the  estimate  is  optimal. 


5. 


THE  CONVERGENCE  RATE  OF  THE  h-p  VERSION.  THE  CASE  OF  THE  SINGULAR 


SOLUTION 

In  Section  4  we  analyzed  the  rate  of  the  h-p  version  when  the  solution 
of  (3-1)  (3-2)  has  the  form  (3.4)  with  u^  =  0.  Now  we  will  analyze  the 
rate  of  convergence  in  the  case  u  =  u^.  For  simplicity  and  without  a  loss 
of  generality  we  will  assume  that  n  =  1  in  (3.4c) 


5.1.  An  approximation  result 

Consider  the  square  R  =  R(h)  defined  in  Section  2.  Let  (r,9) 


denote  the  polar  coordinates  with  the  origin  at  0  (see  Fig.  5.1).  For 
<  >  1  let  S<  be  the  subset  of  R  bounded  by  the  lines  L^:  x2  +  h  * 
<(x1+h)  and  L^:  x1  +  h  =  <(x2+h).  Let  s£  be  the  region  S^  Cl  {(r,0)  | 
r  <  p}  (0  <  p  <  |). 

We  will  consider  the  approximation  of  a  function  u  with  support  in 


s£  for  some  <n  >  <  which  vanishes  on  the  lines  l\.,  L^.  We  will  assume 


that  the  function  u  has  the  form 


(5.1  ) 

u(r,9)  =  ra| 

log  Xo(^)*(0) 

where  5 

and 

Xq  are  sufficiently 

•  smooth  functions  (e.g. 

C* 

functions ) 

such  that 

0  i 

x0  i 

1 .  X0(r)  =  1 

for  0  <  r  <  |,  Xq (r‘)  = 

0 

for  r  >  |-£, 

o  <  p  <  \ 

and 

*(9, 

II 

*> 

CD 

f\J 

II 

O 

where  6-j  ,  02  are  polar 

coordinates  of 

the  lines 

>4 

anad 

■4- 

1  2 

Let  Q  be  the  region  bounded  bv  the  lines  L„  ,  L„  ,  and  x. 

<0  <o  1 

h  h 

~  ~  2’  x2  ~  ~  2  '  We  estimate  the  approximation  error 

|j-uhL  c  pl(R). 


Fig.  5.1.  Scheme  of  R(h),  Q,  s£  . 

IC0 

We  first  map  R  =  R(h)  onto  the  square  R  =  R(1)  by  the  transformation  x 
xi  * 

=  —  or  equivalently  (r,9)  =  (£, 9).  This  maps  Q  into  Q.  Then,  if 

a 

u(r,0)  =  u(r,0)  we  have 

(5.2)  u(r,9)  =  hara|loghr|Y  Xq(^)$(0) 

where  u  =  0  on  the  lines  lJ.  and  L^,  the  maps  of  and  L^.  Since 

y  is  by  assumption  a  positive  integer,  we  have  for  h,  r  <  1 


Y  Y 

(5.3)  u(r,8)  =  T  C(  l)hara|  log  h |  ^ | log  r  |  f  ^\Ar)^(.d)  =  £  uj.. 

2=0  u  2=0 

3y  Theorem  5.1  of  [5]  there  exists  z£  €  Pp+2(R)  such  that  2*  =  0  on  the 
lines  L<  and  and 


I-,  -^1  -  < 

M2  “pJl  ,Q  - 


h^_  |  log  hi  |  |  log  pj  Y  1 

2a 

p 


Hence,  we  see  that 


Zp  =  0  on  and  and 


(5-4)  lu-zpl1  q  <  c  ^  llog  h|i  I108  p|Y  ^ 

,w  p  a  1=0 

ha  V  Y 

<  c  max(|log  h|  ,  I  log  p  J  ). 

P 

By  suitably  changing  the  constant  in  (5.4),  we  see  that  we  may  obtain  a 
*  1  * 

zp  €  Pp(B)  satisfying  (5.4).  By  Lemma  4.2  the  same  estimate  holds  for 
|u-Zp  q  so  that  we  have 

Lemma  5.1.  Let  u  be  given  by  (5.1).  Then  there  exists  Zp  € 
Pp(R)  such  that  zp  =  0  on  the  lines  L^  and  L^  and 

(5.5a)  lu-zpll,Q  <  Cg(h.p,Y)  ^ 

P 

where 

(5.5b)  g(h,p,Y)  =  max(|log  h|Y,  |log  p|Y) 

and  C  is  a  constant  independent  of  p  and  h. 

5.2.  The  rate  of  convergence  of  the  h-p  version 

We  now  return  to  the  problem  of  approximating  the  function  u^  g: 
in  (3.4d).  To  this  end  let 


where 

(5.6a) 


U3 


u3, 1  +  u3 , 2 


*3.1 


=  u 


3  xn(ir) 


-  v  -  v  wvw c-  mt  mywv  »  "T  w  v  —  ’ 


7 ■r».»T.u ■*  v "t  v  *.t ■  r"*  v  v ~y " r»  c -  v  ■ 
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fc -.*  • 

V". 

M 

>*.■-■? 

■•-■■.■ 

i:-“o 


(5.6b) 


u3,2 


u3^ 1  _x0^h^ 


Obviously  2  =  0  in  the  neighborhood  of  the  origin. 

Our  first  goal  is  to  approximate  u^  over  the  set  of  triangles  or 
parallelograms  having  a  vertex  at  the  origin  as  shown  in  Fig.  5.2. 


t.y : 

if.-- . 


m 


rwm 

r.-v  r 

r.%  »*’ 
tv:- 


Fig.  5.2.  Scheme  of  the  mesh  in  the 
neighborhood  of  the  singularity. 


m 


We  will  assume  that  OB^  c  r  and  OB, 


m+1 


c  .  Let 


r  -  U  BiBi+i  •  Then 
i  =0 

Lemma  5.1  yields  the  following  result,  the  proof  of  which  may  be  found  in 


[5]. 


k.V.- 

*  *  « 

L »  *  •  • 

n.%V* 

Kv. 

•  * 

o 

*-* 

X 

o 

& 

5:^ 

exists 

& 

ft* 

(5.7a) 

Lemma  5.2.  Let  u  be  given  by  (5.oa)  with  p  (in  the  definition 


|U-Z  I  <  Cgth.p.Y) 

■  h  <  a)  p 


L  V  %\  S.  /,  ^  ^  ^ 

v'  *.*  *\*  v"  *-'V  %  V' 

Tj-J  >V-( 
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(5.7b)  g(h ,p , Y)  =  max( j log  h  |  Y,  |logp|Y) 


where  C  depends  on  a,  t  but  is  independent  of  p  and  h. 

Let  us  consider  now  the  function  u  =  113  j  given  by  (5.6b).  We 
have  u  =  0  for  r  ^  ph.  Further, 


J  D6u  I  *  C(S)r°H8l  |  log  r|Y 
where  3  =  3^  >  0,  B^  +  B2  =  |Bj  and 


Hence  we  have 


D8u 


Js| 

3  1  1  U 

B1  B2* 
3x1  3x2 


(5.8)  |u|k  Q  <  C(k)|log  h|Y  max(1 ,h“+1  k). 


Denoting  by  u(?  the  finite  element  approximation  of  u,  we  get  by  Theorem 
4.8 

(5.9)  l'‘-“pl,,Q  <  =<*>  Mk.a 

P 

n-k+a 

<  C(k)  p— p  llo«  hl 
P 

with  k  >  1  arbitrary  and  n  =  min(p+1,k).  Let  us  take  k  =  2a  +  1  in 

(5.9) .  Then  n-k+a=n-a-1  =  min(a.p-a)  so  that 

min(a,p-a) 

(5.10)  |’*-u  |  <  C  - p -  I  log  h|  . 

^  *  r\ 


If  p  is  snail  with  respect  to  a,  we  can  select  k  so  that 


C(k)hrrk+Qt/pk~1 


will  be  minimal.  For  example,  with  k  =  2  we  get 


(5.1U  lu"Upll,Q  -  Ch°  | log  h|Y. 

Combining  the  estimates  for  ^  and  2  we  get 

Theorem  5.3-  Let  u  be  given  by  (3.^d).  Then  there  exists 
i>p  €  l/p(£l)  such  that 

.  min(a.p-a) 

(5.12a)  |u-<)>  |  n  ^  Cg(h,p,Y)min(ha,  - 5 - ) 

0  1  m  U  dQL 

P 


(5.12b)  g(h,p,Y)  =  maxCjlog  h|y,  |log  p|y) 

and  C  depends  on  a,  t  but  is  independent  of  p  and  h. 

Remark  1.  When  a  is  an  integer  and  Y  *  0,  the  estimate  (5.12a)  is 
very  pessimistic,  since  the  solution  u  given  by  O.^d)  is  smooth.  When 
a  is  an  integer  and  Y  >  0,  then  the  estimate  (5.12a)  is  a  correct  one. 

Let  us  now  summarize  in  one  theorem  the  error  estimate  for  the  h-p 
version  with  quasiuniform  mesh  and  uniform  p. 

Theorem  5.^.  Let  Q  be  a  polygonal  domain  as  introduced  in  Section 
2.  Suppose  that  u,  the  solution  of  (3»1)-(3.2)  can  be  written  in  the  form 
(3«*0.  Assume  further  if  1  <  k1  i  3/2  that  a  is  a  Lipschitz  domain. 
Assume  that  u^  is  the  finite  element  solution  with  triangular  and 
parallelogram  elements  satisfying  (3.6)  and  the  boundary  condition  on  r1 
defined  by  (3.7).  Then 
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and  C  depends  on  t,  o  in  (3.4),  ft,  k^,  Y^,  but  is  independent 
of  Th,  h,  p,  u, 

Remark  2.  We  formulated  Theorem  5.4  only  in  the  frame  of  Sobolev 
spaces.  By  interpolation  arguments,  it  is  also  possible  to  formulate  the 
theorem  in  the  frame  of  3esov  spaces. 

Remark  3.  We  addressed  only  the  case  of  the  polygonal  domain  and 
elements  which  are  triangles  or  parallelograms.  By  the  standard  mapping 
approach,  the  results  are  also  valid  for  curvilinear  elements. 


6. 


APPLICATIONS 


In  this  section  we  will  study  the  consequences  of  Theor  m  5.4  in 
connection  with  computations. 

First  let  us  mention  that  although  we  discussed  the  h-p  version  in 
connection  with  the  problem  (3.1)  (3.2),  all  conclusions  are  valid  also  for 
the  elasticity  problem.  In  (3.4d)  we  assumed  that  ot^  are  real.  In  the 
case  of  elasticity  problem,  are  in  general  complex  with  Re  a,  >  0. 

The  estimate  (5.13)  is  still  valid  with  =  Re  . 

Our  theory  is  of  asymptotic  character.  Hence  it  is  important  to  see 
the  applicability  of  Theorem  5.4.  in  the  range  of  practical  parameters.  To 
this  end  let  us  consider  the  plane  strain  elasticity  problem  when  ft  is  an 
L-shape  domain  shown  in  Fig.  6.1, 


B  A 


Fig.  6.1.  L-shaped  plane  elastic  body. 

Let  is  assume  that  on  9ft  tractions  are  prescribed,  i.e.  r‘  =  0.  The 
solution  of  this  problem  is  the  displacement  vector  (u-j,U2)  where 


-  id.  A  &  H  a  ^  a  * 


(6.1a) 


—  ra[ ( <-Q(a+1 ) )cos  aO  -  a  cos(a-2)0] 


ui  = 

(6.1b)  u2  =  rQt^ <+Q(a+1  ))sin  a©  -  a  cos(a-2)0] 

where  a  =  0.544483737 

Q  =  0.543075579 


G  is  the  modulus  of  rigidity  and  <  =  3  -  4v  where  v  is  Poisson's  ratio 
which  we  assume  to  be  v  =  0.3.  The  solution  has  a  typical  singularity  at 
0.  The  sides  0A  and  0E  are  traction  free.  Instead  of  the  norm  I*l-|  ^  we 
will  be  interested  in  the  energy  norm  I’lg  which  is  equivalent  to  the 
|,J1  g  norm.  Denoting  W(u),  respectively  W(up),  to  be  the  strain 
energy  of  the  exact,  respectively  the  finite  element  solution,  we  have 


(6.2) 


lu*uolE 


(W(u)-W(u£))l/* 


and  we  define  the  relative  error  in  the  energy  norm  as 


(6.3) 


W(u)-W(u 
^  W(u ) 


h}  \ 

*-]  • 


In  the  next  figures  we  will  present  the  results  of  computations  which  were 
performed  with  a  computer  program  called  PROBE  [20],  [22]  developed  by 
Noetic  Technologies  Corporation,  St  Louis. 

We  will  consider  a  uniform  mesh  with  square  elements  as  shown  in  Fig 

6.2. 


•  -r  ■■-  /.'/.V  V-V- 


A* 

■  W-  /- 


^  ;  ;  •  y-w.  >:v;  .  : 

*j\  * mJ*.  »_r.  \  _v  _•»  J»  _%  a  *j  jI 


Fig.  6.2.  The  scheme  of  the  uniform  mesh. 

The  solution  u  €  H1+ot-e(n),  e  >  0  arbitrary. 

Theorem  5.4  gives  for  p  >  1  the  estimate: 

min(a.p-a) 

(6.4)  |u-u  |E  £  C  min[ha,  - jjjj - ] 

p  p^a 

where  C  depends  on  a  but  is  independent  of  h  and  p.  Fig.  6.3  shows 
the  relative  error  in  the  energy  norm  |e|r  Q  (for  different  degrees  p) 

U  |  11 

in  dependence  on  h.  We  also  show  the  slope  ha  in  the  figure.  We  see 
that  with  respect  to  h  the  error  is  in  the  asymptotic  range  also  for 


moderate  p  and  h. 


rig.  6.3.  The  relative  error  in  the  energy  norm  in  dependence  on  h. 


Fig.  6.4  shows  the  error  in  dependence  on  p  and  different  h.  Because  of 
the  size  of  computations,  only  in  the  case  h  =  \  is  the  error  given  for  p 
>  4.  (For  p  =  4  and  h  -  1/10,  the  number  of  degrees  of  freedom  N  = 

5119).  Estimate  6.4  gives  the  rate  P~^a  which  appears  only  for  p  >  3. 

2  2 

For  large  p  and  small  h  we  have  N  =  p  /h  and  hence 


(6.5) 


a 


(6.5)  shows  that  if  the  measure  of  computational  work  is  N,  then  the  use 


of  higher  p  is  preferable. 


IN  ENI 
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Fig.  6.5  shows  the  dependence  of  the  relative  error  in  the  energy  norm 
on  the  number  of  degrees  of  freedom  N  for  various  p.  In 
addition,  the  performance  of  the  p-version  for  h  =  V2  is  shown  in  the 
figure.  We  see  that  p  =  2  is  more  effective  than  p  =  3,  and 
asymptotically  for  p  ■*  ”,  the  higher  p  are  more  effective  as  follows 
from  (6.4).  The  p-version  has  a  rate  which  is  twice  that  of  the  h-version 
(see  also  [5]) . 

We  addressed  in  this  paper  only  the  case  of  the  quasiuniform  mesh.  If 
the  mesh  is  strongly  refined,  then  its  performace  is  different.  Fig.  6.6 
shows  the  strongly  refined  mesh  with  n  layers,  (n*2).  The  mesh  is  a 
geometric  one  with  the  ratio  0.15.  The  ratio  0.15  leads  to  nearly  optimal 
convergence.  See  C 1 3 1 »  [14]. 

Fig.  6.7  compares  the  performance  of  the  h,  p  version  for  the  uniform 
and  strongly  refined  mesh  for  our  example.  The  performance  of  the  p-version 
on  strongly  refined  meshes  is  in  practice  very  similar  to  the  general  h-p 
version,  leading  to  an  exponential  rate  of  convergence.  We  see  that  the  p- 
version  performance  depends  very  strongly  on  the  mesh. 

For  more  about  the  comparison  between  h,  p  and  h-p  version  we  refer 
to  [3  ]  - 


NUMBER  OF  DEGREES  OF  FREEDOM 


<  o 


d  mesh  with  n  =  2  layers. 


Ml 


400  8001000  2000  4000  0000 


NUMBER  OF  DEGREES  OF  FREEDOM 


rig.  6.7.  The  error  in  the  energy  norm  in  dependence  on  N  for  various 
meshes . 


7.  APPENDIX 

Theorems  7.4  and  7.5  proven  in  this  section  are  slightly  generalized 
forms  of  Lemma  4.7  and  are  of  interest  by  themselves. 

Let  us  consider  the  equilateral  triangle  T  =  ABC  as  shown  in  Fig. 
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(7.1  ) 


F1Cf](x,y) 


x+  JL- 
/3  f  ^3 

g-  J  f(t)dt. 

x-  i- 
/3 


The  value  of  F-j  at  a  point  P  €  T  depends  only  on  the  values  f  along 

the  segment  Q.Q-.  Qi  =  (x-  ,0),  Q~  =  (x+  — 0).  We  prove  now  the 

12  /3  2  /3 

following  lemma. 


7.1.  Let  f  €  P  (Y-j)  and  F^f-*(x,y)  be  defined  by  (7.1) 


a) 

pff](x,y)  €  Pj(T) 

b) 

Fff](x,0)  =  f ( x ) 

c) 

Clfl-4.Y, 

di> 

lFff]l  A  i 

k ,  Y^ 

c|f|  A 

k.Y, 

0  * 

k  *  1 

d2) 

|rff]I  B  i 

k ,  Y^ 

clfl  B 

k,Y 

o  < 

k  S  1 

d^) 

lF1Cf]|  c  < 

Clflo,Y, 

0  1 

k  £  1 

d4) 

lF1Cf]I  r  < 

k,Y" 

Clfl0.Y. 

1 

0  1 

k  £  1 

where  the  constant  C  is  independent  of  p  and  f. 
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Proof.  It  is  immediate  that  (7.2b)  holds.  Let  f  =  xn  with  0  <  n  £ 


p  integer.  Then 


F(x,y)  = 


x+ 

0  Av 

v-  L_ 


-JL _ [(x  ♦  L_)n+1  -  (x  -  *-)n+1] 

2y(n+1  )  UX  ^  yj  J 


2yfn^TT  *  (x-  ‘  TSW  V^>  5 


Hence  (7.2a)  holds. 

To  prove  (7.2c)  we  first  extend  f  to  a  function  defined  on  the 
entire  x-axis  R  so  that  (see  [18] ) 


(7.3) 


S  clfl-4 


where  we  have  used  the  same  notation  f  to  denote  the  extended  function  as 
well.  Then  by  (7.1)  F^x.y)  is  well  defined  on  the  entire  half  plane  J2  = 
{(x,y)  |  y  >  0}.  For  (x,y)  €  fl  we  have 


(7.4) 


F.,(xfy)  =  j  f (t)H(x-t,y)dt  =  (f  *  H(*,y))(x) 


where 


(7.5) 


H(x,y)  = 


-  <  x  < 

/3  "  "  /3 


0  otherwise. 


Let  g(C)  represent  the  Fourier  transform  of  the  function  g(x)  in  the 
x  direction.  Then  by  (7.4) 


(7.6)  FjU.y)  =  f(OHU.y) 

where 


(7.7) 


HU.y) 


J_/3  if 

/2?  2y  -y//3 


-i  £x 

e  dx  = 


1  sin(  C.y//3) 
/2ir  £y//3 


Let  ft  =  {(5,y)  |  y  >  0}  and  calculate  the  H^Q)  norm  of  F^x.y).  By 
Parseval’s  equality,  we  have  using  (7.6) 


lF!  I2!  =  IM2,  _  =  \\  |fU)|2U  H(5,y)|2d5  dy 

H  (0)  H  (Q)  - 

*  JJ  |f(0|2  ||^r  HU.y)|2dj:  dy  ♦  JJ  |f(U|2  |HU,y)|2d5  dy. 
a  a 


Now  letting  z  =  y£//3  we  get,  by  (7.7), 


03 


|  |HU.y)|2dy 
0 


00 


sin  z  dz 

~T?T  * 


C 

ITT  ‘ 


Hence 

(7.9)  JJ  [f(C)|2UH(5,y)j2dC  dy 
a 

ao 

<  Cj  |5||f(5)|2d5  <  C|f|2  <  c|f|2y  . 

-a,  '*»  "  /2»  '  i 

Also 


.  AM  lt«V  A  ,^kV  A  -  -V  ■.’V.% 
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!_fiu,y)  .  ^_:=2p 


sin  2- 


which  is  bounded  at  2=0.  Hence 


so  that 


J  ih  S(5,y)|2dy  <  C|5|, 


(7.10) 


SI  |fU)|2  if-  H(C.y)  |2d^  dy 


<  C  J  |c||f(5)|2dC  <|f|2  y 


The  third  term  can  be  bounded  analogously.  Using  (7. 8)- (7-1 0) ,  (7.2c) 
follows.  Inequalities  (7.2d^)»  (7.2djj)  follow  immediately  for  k  =  0, 
k  =  1  and  hence  by  an  interpolation  argument  (see  [8])  they  hold  for  all 
0  £  k  <  1  . 

We  prove  now  ( 7 -  2d -| )  -  Let  the  variable  x  be  used  to  represent  both 
the  distance  from  A  along  Y1  and  the  distance  from  A  along  >2* 
Denoting 


(7.11  ) 


1  x 

G(x )  =  -  f  f (t)dt 

Y  J 


it  is  readily  seen  that 


(7.12) 


|F,Cf]| 


|G(x)|, 


I  =  (0,l/2). 


Using  (9.9.1)  of  [15  ],  p.  244  we  get 


Further,  integrating  (7.11)  by  parts  we  have 


(7.14) 
and  hence 


G(x)  =  f(x)  -  f (0) 


-iJ 


tf '(t)dt 


G'(x)  =  f’(x)  +  j  tf '  (t)dt  -  f’(x) 

x  0 


.  x  .x 

-=■  I  (x-t)f ' (t)dt  ♦  —  l  f * (t)dt. 

/  ^  A 


x2  0 


Using  9.9.5  of  [15],  p.  245  with  r  =  2  we  get 


1^  I  <x-t)f'(t)dt|0il  <  C|f|0I 


and  by  9.9.1  of  [15],  p.  244  we  get 


\h  I  f’(t)dt|0iI  <  C|r| 


'o,r 


Hence 

(7.15^  3gt(x)30,i  <  Clf10,I  • 


Conbining  (7.13)  and  (7.15)  we  get  (7.2.d^ )  for  k  =  0  and  k  =  1 
hence  by  the  interpolation  argument  ( 7 -  2d i  )  holds  for  all  0  <,  k  <,  1 
inequality  (7.2d2)  is  essentially  the  same  as  (7.2d1)  an(j  Lemma  7.1 


completely  proven. 
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Let  now  f  =  f^  €  P  (Y^),  i  =  1,2,3.  Then  we  denote  by 

[f,] 

F^  (x,y)  the  polynomial  extension  of  f<  into  T,  defined  for  i 
by  (7.1)  and  for  i  =  2,3  by  (7.1)  after  properly  rotating  the 
coordinates.  Obviously  Lemma  (7.1)  is  applicable  for  i  =  1,2,3  when 
properly  interpreted  through  the  rotation  of  the  coordinates. 

We  now  prove 


Lenma  7.2.  Let  T  be  the  triangle  as  in  Fig.  7.1  and  f  be  a 

continuous  function  on  3T,  such  that  f ^  =  f|Y  €  Pp(Y^),  *  =  1*2,3 

where  by  f|Y  we  denote  the  restriction  of  f  on  Y. .  Then  there 
i 

exists  €  P  ( Y . ) ,  i  =  1,2  such  that 
1  ^  1 


[♦,]  [*,]  „ 

(7.16)  a)  U  =  F}  1  +  F2  d  €  P  (T) 


b)  U  =  ft  on  Y,,  i  =  1  ,2 


G)  lUll,T  -  Ct-Ifll‘/2,Y1  +  Hf2ll4,Y23 

V  i$ilk,Y.  *  C[.^  lfjlk,Y.3' 

d2>  JM  B  -  C[lfll  B+  *  lfjl0,Y 

k,Y1  k,Y“  J-1  J  j 

2 

d3)  3^1  C  *  CC  lr2*  C  +  1  Ir i  l0  y  3. 

*  k, Y^  k , Y^  J-l  J  U*  j 


i  -  1 ,2,  0  *  k  £  1 


0  £  k  <  1 


0  *  k  i  1 


where  C  is  a  constant  independent  of  p  and  f. 


Proof.  Let  €  P  (Y^).  Then  as  in  Lemma  1  we  define 


r  ./../• 


'.y.  v 

.•  /  /  *  *  1  •  ’  •  ’  ••  »  * 


1  *  •„  1  >‘v\*  -‘..v  V,y;  >V.  \% 

■  V  -  \  "-’v  V  V  V.  *u\ »JV  k.*V\\L.V  A 


(7.17) 


i  =  1.2 


1  A 

GjCx)  =  -  J  *. (t)dt, 
x  0  1 

Condition  (7.16b)  will  be  satisfied  if 

(7.18a)  (x)  +  G2(x)  =  (x)  +  —  J  4>_(t)dt  =  f  (x) 

x  o  ^ 

(7.18b)  *2(x)  +  G^x)  =  *2(x)  +  7]  ♦1(t)dt  =  f 2 ( x ) 

hold  for  all  x  €  I  =  (0,1).  Since  fi  c  Pp( I)  it  is  easy  to  see 

that  ^  €  Pp(I)  satisfying  (7.18)  exist.  Due  to  the  assumption  about 

continuity  of  f  on  3T  we  have  f 1 ( 0 )  =  f2(0)  =  C.  are  uniquely 

determined  up  to  a  constant  K  with  (0)  =  K,  $2(0)  =  C  -  K. 

We  now  define 


(7.19) 

(x) 

*  ♦1  (x)  +  lf2(x) ,  ^(x) 

*  ^(x) 

-  $2(x) 

h1  (x) 

=  f1 (x)  +  f2(x) ,  h2(x) 

=  fj(x) 

-  f2(x) 

30  that 

(7.18) 

yields 

(7.20a) 

1  x 

'J'l  (x)  +  7  f  (t)dt  = 

0 

h1  (x) 

ip9(x)  -  7  ?  1 Mt)dt  =  h0(x). 

6 


(7.20b) 
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Here  '^(x)  is  unique,  ^  (0)  =  C,  while  i|>2(x)  is  unique  up  to  the 
constant  K  such  that  <|>2(0)  =  2K  -C. 

We  first  analyze  (7.20a).  By  differentiation  we  obtain 


~  1  <l>i  (t)dt  +  ^ 


=  h,  . 


Using  (7.20a)  we  get 


(7.21  ) 


x 


h*  +  — 
1  x 


The  homogeneous  solution  of  (7.21)  is  — .  A  particular  solution  can  be 

x 

found  by  using  the  method  of  variation  of  constants.  Hence,  substituting 

(x)  =  into  (7.21)  we  get 

x 


T'(x)  =  h{x2  ♦  hlX 


from  which 


*-|  (x) 


=  —p  /  t2h1T  (t)dt  +  \  th^Odt. 


Integrating  by  parts  we  get 


(7.22) 


^  (x) 


=  h1 (x)  -  J  th1 (t] 


the  unique  solution  of  (7.20a) 


We  show  now  that 


hilk.i  ^  cMk>r 


0  <  k  <  1. 


,» V  / 
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where 


A  A  « 

j  th1(t)dt  =  -  j  (x-t)h1(t)dt  +  x  }  h1(t)dt. 

0  0  0 


F(x )  G(x ) 

2  =  2 
X  X 


-  Q(x ) 


G(x)  =  J  (x-t)h1(t)dt 
0 


Q(x )  =  j  h1  (t)dt, 

0 


Using  (9.9.4)  of  05],  d.  245  with  r  =  2  and  (9.9.1)  of  [15],  p.  244  we 


obtain 


|F(x)/xa|0>I  <  |G(x)/x2|0J.  |q|0jI  <  C|h,l0J 


which  yields  (7.23)  for  k  =  0.  Next,  differentiating  (7.22)  we  get 


(7.24)  ^ 


P  yr  h1  1  xr  2 

h{  +  ^  j  th1  (t)df  -  —  =  hf  -  -j  JVh*(t)dt. 

x3  0  x^  0 


X  X 

F(x )  =  /  t2h'(t)dt  =  J  (x-t)2h’(t)dt 

0  1  0 


-x2  f  hj(t)dt  2x  ]  th’(t)dt. 


We  have  then 
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where 


F(x)  =  G(x) 
x3  x3 


-  Q(x )  +  R(x) 


A 

G(x)  =  { 


(x-t)  h.j(t)dt 


Q(x) 


x  i  1 


,'(t)  dt 


This  gives 


R(x)  =  %  J  th'(t)dt. 

x  0 


JF(x)x3  |0j  j  <  |G(x)x3|0tI  ♦  5Q(x)||0>i  ♦  |R(x)I0J, 


The  first  two  terms  can  be  bounded  once  more  by  |hj |q  j  using  (9.9.4)  of 
[15],  p.  245  and  (9.9.1),  p.  244.  Moreover, 

R(x)  =  ^  [-/  (x-t)hj(t)dt  +  x  /  h' (t)dt] 

*  0  0 

so  that  |R|0  j  can  also  be  bounded  by  l0  j*  This  yields  (7.23)  for  k 
=  1.  By  the  interpolation  argument  (see  [8])  we  get  immediately  (7.23).  Let 
us  consider  now  (7.20b).  Differentiating  it  and  using  once  more  (7.20b)  we 
get 

*  «  ^2 

(7.25)  ^2  =  ^2  +  ~  ’ 

Integrating  we  get 

1  h  (t) 

(7.26)  4*2 ( x )  *  h2(x)  -  j  — - —  dt. 

x 
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(7-26)  is  that  solution  of  (7.20b)  with  ^7(1)  =  h0(1). 

Once  more  we  wish  to  show 

(7-27)  IMk.I  -  C  1^2  Ik ,T  ’  0  *  k  *  1 

Using  (7.26)  and  (9.9.9)  from  [15),  p.  245  with  a  =  0  we  get 

l^lo.I  -  u  ®'n2  Iq , T  * 

Since  h2(0)  =  f^O)  -  f 2(0)  =  0,  (7.25)  yields 


I 


1  A 

h2  +  -  J  h’(t)dt 


and  by  (9.9.1)  of  [15],  p.  244  we  get 

l^ll  ,  J  -  ,  7  • 

An  interpolation  argument  leads  immediately  to  (7.27).  Hence  we  have 
constructed  solutions  of  (7.20a,b)  such  that  (7.23)  and  (7.27)  hold. 
Coming  back  to  (7.19),  using  k  =  \  we  see  that  for  i  =  1 ,2 


i*t  l* 


<  C[|f  |y 

1  '2»  1  -j 


I r  p  I  y  v  ^ 

c.  /2,  1 2 


and  applying  Lemma  7.1  we  get  immediately  (7.16c)  and  also  (7.1 6d1 ) . 
Returning  to  (7.20)  we  see  that  with  I*  =  (l4,1) 


which  immediately  leads  to  (7.l6d2).  (7.1 6d 3 ) - 
The  following  lemma  is  taken  from  [6] 

Lemma  7.3.  Let  T  be  the  triangle  as  before,  f  be  continuous  on 
3T,  f2  =  f^  =  0  and  f1  €  PpC^).  Then  there  exists  a  polynomial 
v  6  Pp(T)  such  that 

lvl1)T  5  clfili,Yi 

v  =  f1  on  Y1 
v  =  0  on  Y2,Y3 

where  C  is  a  constant  independent  of  f  and  p. 

Theorem  7.4.  Let  T  be  the  equilateral  triangle  shown  in  Fig.  7.1 

and  f  be  a  continuous  function  on  3T,  such  that  fj  =  f|Y  €  Pp(Y^), 

i  =  1,2,3.  Then  there  exists  U  €  P JL ( T)  such  that  U  =  f  on  3T  and 

P 

H  T  <  C[  I  |f  ly  y  ] 

’  1  i  =  1  4,Ti 

where  the  constant  C  is  independent  of  p  and  f. 

Proof .  Using  Lemma  7.2  we  see  that  without  loss  of  generality  we  can 
assume  that  f2  =  f^  =0- 

Let  f1  4  0,  f2  =  0.  By  Lemma  7.2  we  construct  $1  ,  $2  and  U  = 

'  ■*,  ]  [<£>?] 

F1  1  *  F2  *  .  Then  U  €  P  (T),  U  =  fj  on  Yj_ ,  i  -  1,2  and 

!,Jll  .7  S 


(7.28) 


Denote  by  g^  the  trace  of  U  on  Y^.  Then  we  have  83(B)  =  83(C)  =  0 


(7.29) 


^3^4,  y3  -  c^fil‘4,Y1 


by  applying  Lemmas  7.2  and  7.1 


Because  of  (7.1  6do )  r  <  C |f ,  |X/  and  hence  using  Lemma  7.1 

3  —  4  v/L  »  'it  '  1 


we  have  also 


(7.30) 


b',  YC  i  C«f1  1, 

1  t  ' 


Let  now  analogously  as  before 


“1  •  h  ♦  ?i 


so  that 


U1 


=  83  on  y3, 


=  0  on  y. 


(  7  -31  1 


IU1 1 1 , t  -  c^s3^4,y  -  clfi  •14,Y1  • 


Denote  by  g^1  J  the  trace  of  on  Y2.  Then  g£IJ(A)  =  g?' J(C)  =  0. 

Because  of  (7.30),  applying  Lemma  7.1  and  Lemma  7.2  analogously  as  before  we 
conclude  that 


[1  .  ,[13/ 


3llfyC  +  *83*,4,y3]  -  c!fi  J'4,y1 


Mow  applying  Lemma  7.3  there  is  U2  e  P^CT)  such  that 


Let  now 


U2  =  S2  on  y2’  U2  =  0  on  Y1  '  Y3’ 

V  =  U  -  u1  +  u2. 

Then  it  is  easy  to  see  that  V  €  Pp(T),  V  =  f1  on  Y1 ,  V  =  0  on  Y2, 
and  because  of  (7.28),  (7.31)  and  (7.32)  we  get 

lvll>T  <  clf1li/2Yi 

which  concludes  the  proof  of  Theorem  7.4. 

Let  S  =  (x.yj  |x|  <  1,  |y|  <  1)  be  a  square  and  Y^  its  sides  as 

shown  in  rig.  7.2 


Fig.  7.2.  The  scheme  of  the  square. 


Theorem  7.5  Let  S  be  the  square  shown  in  Fig.  7.2  and  f  be  a 
continuous  function  on  3S,  such  that  f^  =  f|y  €  PpC"^),  i  =  1  .-..,4. 
Then  there  exists  U  €  Pp(S)  such  that  U  =  f  on  3S  and 


where  the  constant  C  is  independent  of  p  and  f. 

Proof.  Let  T  be  triangle  shown  in  Fig.  7.1  and  Q  *  (5,n  | 
(£,n)  €  T,  n  <  ^^}be  the  trapezoid  shown  in  Fig.  7.3. 


Fig.  7.3.  Scheme  of  the  trapezoid. 


The  mapping 

(7.33)  5  =  \  +  If  ('y  +  3°  n  =  (1  +y)  * 

maps  3  onto  Q.  The  mapping  is  obviously  one-to-one  and  the  Jacobian  and 


it3  inverse  are  bounded 


Let  us  first  prove  the  theorem  in  the  case  that  f^  =0,  i  *  2,3.^. 
Denote  f1  U)  =  f1(2^-1),  0  <  £  <  1. 

Obviously  f-|(A)  =  f.j(B)  =  0  and 


Let  0  €  Pp(T) 
Theorem  7-4, 


l?ll  - 

1  1/  Y 


y  Y 

'2*  'i 


such  that  U  =  f 


1 


U(5,n)  exists  and 


£  C  l^i  |y  -y  • 
on  55  and  0=0 


on  5£  and  5C. 


By 


H,t  <  c|f,|  <  c|f,  |,4>  . 

1 


Because  0  €  Pp(T)  we  have 


U(5,n)  = 


XJ 


I  ak  Crt 

0<k+j<p 


I  "  rl  C-y*|)>K(<i-y)  S£I_)J 


o<k+j<p 


=  U(x,y)  €  Pp(S) 


and 


l"l,,s  i  clfl'4,r,- 


Because  f2  =  ?3  =  0  we  have  U(±1 ,y)  =  0,  U(x,-1)  =  f1  and  using  Lemma 

7.1,  Lemma  7.2  we  conclude  by  similar  arguments  as  used  in  the  proof  of 


Theorem  7.1*  that 
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Of  course  U(x,1)  €  Pp(Y3)  and  U(±1,1)  =  0.  Hence  with 

V  =  V2  U(x,1)(y+1 ) 

we  see  that 

|V|,,S  <  C|U(x,t)|tjY3  <  |f,|.4>Y] 

and  V(x,1)  =  U(x,1).  Hence  W  *  1)  -  V  (  P^(S),  W  =  f  on  3S  and 

l«li  .s  s  clfi 


The  theorem  is  therefore  proven  in  the  case  that  f  =  0  on  three  sides  of 
S  and  hence  it  holds  also  if  f  is  general  but  f  =  0  at  the  vertices 
ABCD. 

It  remains  to  prove  that  in  the  general  case  there  exi3t  4>  €  Pp(S) 
such  that  4  has  the  same  traces  at  ABCD  as  f  and 


(7.34) 


Mi.s  -  C^1  lfily2,Y.^ 


Cf,]  [f,] 

To  this  end  we  define  F1  (5,n)  by  (7.1)  and  define  F1  (x,y)  by 

[f1] 

inserting  (7.33)  for  (£,n).  Then  |F.  (x,+1)|.  y  <  cjf.  |n  Y  and  hence 

1  1 »  '3  1  u> 

Cf,] 

analogously  as  above  we  can  change  F^  (x,y)  to  F1  so  that  F(x,+1)  =  0 

and  JF^  ]1  s  <  C | f 1  |y  y  .  Changing  the  role  of  Y1  and  Y^  we  can 
analogously  construct  F^  $  P^(s)  so  that 


|f3|1>s  <  c|f3|y^Y  ,  F3(x,1)  *  f3,  F3(x,-1)  =  0. 
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